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Abstract: A simple model, based on two parallel and independent NN —>■ dn processes, 
has recently been proposed for two-pion production in the dd — >■ aX reaction. It reproduces 
all observed features, including the sharp peak structure in momentum distributions (the 
ABC effect) and the strong oscillations in the deuteron vector and tensor analyzing powers. 
This model is now extended to describe also four-pion production with the same basic 
mechanism, but with two np dTnr processes as input. The calculations of the high 
missing mass spectra are within about 30% of the experimental data for beam energies in 
the range 1.9 < < 2.4 GeV. 
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It has long been known that the momentum distributions of the np dX 12 E] , 
pd — ^HeX in |3], and dd aX ^ [7j reactions at intermediate energies share some 
common and spectacular features. Characteristic are the sharp peaks, located slightly 
above the two-pion threshold, in association with a broader bump around the maximal 
missing mass. This behaviour, known as the ABC effect, is observed only in isospin Jx = 
channels. The phenomenon is believed to originate from a largely kinematic enhancement 
in the production of two pions, but all previous attempts to give a quantitative description 
have so far had only limited success. In particular, the theoretical models for the np —>■ 
dX [HI ini Hm mj and pd ^B.eX jT21 HHj reactions have problems in reproducing the 
angular dependences [SI [12] • A recently published model for dd aX [13] has, however, 
had remarkable success in reproducing data. It is based on two almost free NN dir 
processes happening in parallel, where the two final deuterons merge to form the a particle. 
Since the deuteron is loosely bound, these processes are almost free and their cm. systems 
coincide with the overall cm. frame. Because of the strong forward-backward peaking 
in the two dn systems, the small missing mass (parallel pions) and maximal missing mass 
(back-to-back pions) configurations are both enhanced. In the latter case, the large relative 
momentum between the deuterons in the a particle does cause a suppression relative to 
the low TTTT mass case. The ABC peaks and the central bump are hence simultaneously 
explained. The model, depicted by the Feynman diagram in Fig. 1, is able to reproduce 
quantitatively all the observed features of dd — > aX in the energy range 0.8 < Td < 
1.9 GeV, including angular and energy distributions [13] and deuteron vector and tensor 
analyzing powers [3 [^ . The purpose of this work is to use the same type of mechanism 
to estimate four-pion production with two parallel np dinT processes as input. 

At higher energies {Td > 1.9 GeV), where production of heavier mesons and multiple 
(> 2) pions begins to be important, the central bump in 27r-production at large tttt effective 
masses is strongly suppressed. The large relative momentum between the two pions requires 
the deuterons to be far out in the tail of the a-particle wave function in momentum space. 
The angular dependences of the subprocesses are here of minor importance — it is the 
dynamics of the dd : a vertex that dominates the behavior. The situation will be very 
similar when heavier mesons are produced by the same mechanism, the effect being more 
pronounced the larger the maximal missing mass is compared to the meson mass. At 
high enough energy the formation of the a particle will force the two mesons to be close 
together in momentum space and suppress other configurations. For the production of 
pairs of mesons with definite masses we then expect sharp peaks close to the two-meson 
threshold and a rapid decrease towards maximal missing mass. On the other hand, when 
the meson has a broad width, like the controversial a meson, the peaks corresponding 
to different masses should be added together with appropriate weights. In the present 
analysis we will use the a meson as a convenient tool to describe and handle the Ix = 
pion pair produced in np —>■ dX. By considering two parallel a productions, the four- 
pion contribution could be estimated in the same framework as the two-pion one. The 
missing mass distributions of np dX is then the weight functions needed as input to 
this dd aoa model. 

Since the mathematical details of the model have been published ^3], only the main 
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features are recapitulated here. The dashed hues in the Feynman diagram (Fig. 1) are now 
considered to be a mesons with the (possibly unequal) masses mi and Instead of the 
complicated and many AW dir amplitudes used in our previous work, we here assume 
that the np ad reaction can be parametrized by one single s-wave amplitude 



Vp, (1) 



where rix is a nonrelativistic two-component spinor, 77^ = —ia2f]* its charge conjugate, 
the deuteron polarization vector, and A = A{m„) the a meson production amphtude. 
After summing over internal spins the spin kernel /C is, with this normalization, 

IC = ^^A\ (2) 
2x/3 ^ ^ 

where are the polarization vectors of the initial deuterons. The summed square is then 
J2 = (1/4)1^1'^, but an extra factor of four is needed because there are two ways of 
distributing the four nucleons between the two ad vertices. 

The mass dependence of the amplitude A was obtained from experimental data for 
np — dX jll 0] by transforming the laboratory momentum distribution to a missing mass 
distribution in the cm. frame, where it should be symmetric. The data were taken with a 
spectrum of neutron momenta and their quality suggests that one should use the forward 
ABC peak (low k]^^) at the highest energy [3|, while the backward ABC is the better choice 
at 991 MeV p. By integrating this over a small mass range Am^ the differential cross 
section for ad could be estimated: 



d'^a 

Anix. (3) 
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In the normalization we are using, the relation between the cross section and the amplitude, 
and in particular its mass dependence, is given by 



' ' Am,, (4) 



16(27r)2 So-d p dm. 



where m is the nucleon mass and p and k are the nucleon and deuteron cm. momenta. The 
squared amplitude d\A\'^ / drrix, extracted from Eqs. (jHI) and (jH), was fitted to a third degree 
polynomial in mx for each beam energy. At Td = 1163 MeV an angular distribution has 
been measured jH] and the transformation was done for each of the angles. The variation 
of the amplitude with the a mass is roughly independent of the laboratory angle, which 
supports our s-wave assumption, and the polynomial was fitted to all of the 1163 MeV data 
simultaneously. Since none of the energies where the amplitude functions were extracted 
coincides with the energies required for our dd estimate, a linear interpolation in Ei = ^/s^ 
was applied between the fitted curves. The masses used in the fitted polynomials were 
estimated by letting them be in the same proportion to the mass range at the fitted 
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energies as the actual mass is to its allowed mass range. The amplitudes were then finally 
evaluated for each individual m^, by multiplying the fitted and interpolated values by the 
increment Arrix- 

The matrix element Aiaaa is written as the product A4 = —i{ma/vd)]CyV, where W 
is the dimensionless form factor defined in Ref. jlSj. In this factorization approximation 
the amplitudes are calculated without considering the Fermi momenta. Nevertheless it 
is possible to get a crude estimate of their infiuence by assuming that the distribution 
of energy between the two subprocesses follows a simple Gaussian, exp[— (Aii^j)^/(2cr|;)]. 
Here AEi = Ei — is the deviation in energy from equal sharing and aE ~ 25 MeV 

is calculated from a Gaussian deuteron wave function. The calculations are very stable 
against variations of this distribution. Only for small assumed widths < 5 MeV) could 
significant changes be observed. 

The calculation of the form factor W follows the lines of ^H] , apart from obvious mass 
and kinematical changes. The integration over meson angles was performed directly since 
there are no angular or momentum dependences in the amplitudes. For each choice of total 
missing mass M^, the momentum distributions were calculated for all possible combination 
of the two 0" masses. The results were then accumulated to give the full four-pion spectrum. 
With this procedure, all the allowed four-body pion phase space was completely accounted 
for. 

Our predictions of the four-pion missing mass distributions in the cm. frame can be 
seen in Fig. 2 together with data points taken from . Note that the scale factors applied 
to the calculations are not directly comparable to those in [13j; the latter values should 
be multiplied by the relativistic factor 1/7^ = (m^/i^rf)^, which was overlooked in the 
programs (but not in the formulae) of that reference. 

A plausible reason for the underestimate of the cross section is the neglect of contribu- 
tions from pp —>■ ppiTTT. The cross section for this process is at ~ 1.2 GeV about half that 
for np dmr but, lacking information about angular distributions and the nnpp : a 
wave function, it is impossible to describe its influence. Another effect will occur at low 
missing masses {Mx ~ 4m7r), where the statistics of 4 equal bosons should increase the 
counting rate. In the semi-classical model we are employing — and with no knowledge of 
the np dim amplitudes — this effect cannot be estimated. Distortion due to rescattering 
of the pions might decrease the cross section. 

Superposed upon all three of the in spectra there is an additional structure which, 
since it peaks at Mx ~ 780 MeV/c^ independent of beam energy, might be attributed 
to the u meson. This interpretation is however not unambiguous since the widths of the 
peaks are much larger than the stated experimental resolution. The experimentalists noted 
this difference and remarked that the resolution was calculated from known properties of 
the beam and spectrometer but that they were never able to test it directly because of 
poor statistics for two-body reactions in dd aX ^Sj. The upper curves in Fig. 2 result 
from single Gaussians fitted to these points which, after integration give the cm. cross 
sections da/dQ{^uj^) = 9, 6, and 2 nb/sr for 1.9, 2.2, and 2.4 GeV. Because of the threshold 
An probable enhancement mentioned in the previous paragraph, these values are probably 
overestimates. 
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The assumption that the four pions are produced via two a's has the consequence that 

the relative weights for charged and neutral pions should be A^(+ H ) : iV(0 H — ) : 

A^(OOOO) = 4:4:1, which would be a test of the model, since other possible production 
channels have different weights. In the case of pp production, where 47r° is forbidden, 

A^(+ H ) : A^(0 H — ) =2:1. This contribution can be estimated from pp — >■ (iyr^TT"'". 

The cross section for this process is smaller than that for np — > c/tttt by roughly a factor of 
10 Pl] and, in addition, there is only one diagram in this case, giving an extra factor 1/4. 
In total this will give a contribution of less than 1% of the aa cross section. 

Despite the crudeness and simplistic nature of the assumptions and approximations 
made in this model, the calculations are in pleasing agreement with data. Taken together 
with the results of JHl, it seems likely that most of the multi-meson production in the 
dd — ^ aX reaction in the 0.8 < < 2.4 GeV range could be attributed to the same 
general double-interaction mechanism. Further support for this suggestion could be ob- 
tained from the estimate of 77 production via two np — ^ (i7 processes. This is shown to 
be of considerable importance as background to the charge-symmetry-breaking dd —>■ air^ 
reaction 

1 am grateful to Colin Wilkin for suggesting this approach and would like to thank him 
and Goran Faldt for stimulating discussions of the results. 
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Figure 1: The Feynman diagram for the dd — >• anrr and dd aaa reactions. 
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Figure 2: The missing mass distributions for dd — > aX at 6'^ ^ 0°. The data 5j are 
compared with the aa (lower sohd hne) model, the latter curves multiplied by scale factors 
1.2, 1.0, and 1.3 for 1.9, 2.2, and 2.4 GeV. The upper solid line is the fit to the 'cj' peak. 
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